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Quantum Chromodynamics is an established part of the Standard Model and an essential part of the 
toolkit for searching for new physics at high-energy colliders. I present a status report on the theory 
of QCD and review some of the important developments in the past year. 



1 Introduction 

Quantum Chromodynamics (QCD), the 
gauge field theory that describes the inter- 
actions of coloured quarks and gluons, is one 
of the components of the SU(3) x SU(2) x U(l) 
Standard Model. At short distances, equiv- 
alently high energies, the effective coupling 
is small and the theory can be studied us- 
ing perturbative techniques. Enormous effort 
has been made over the past three decades to 
calculate higher-order pQCD corrections to 
phcnomcnologically relevant quantities and, 
as we shall see in this review, the fron- 
tier is now at next-to-next-to-leading order 
(NNLO). Non-perturbative QCD contribu- 
tions are also important. They control, for 
example, the transitions between hadrons 
and quarks and gluons, both in the initial 
and final states, and thus have a key role to 
play in dictating the overall 'shape' of a high- 
energy collider event. Again, there has been 
much progress in understanding and mod- 
elling these non-perturbative effects in re- 
cent years, including approaches based on lat- 
tice calculations, Regge theory, Skyrme and 
large- iV c models etc., to the extent that in 
many areas of application QCD phenomenol- 
ogy is now a high-precision science. How- 
ever there are still gaps in our understand- 
ing, in particular there are processes for 
which our knowledge is still only at the semi- 
quantitative level and much more work needs 
to be done. Examples are semi-hard, exclu- 
sive and soft processes at colliders. This is 



another important frontier for QCD. 

'QCD as a high-precision tool at collid- 
ers' is the theme of this talk. I will begin 
by reviewing the status of pQCD calculations 
and as measurements, and then highlighting 
some significant calculational developments 
in the past year. This includes the emergence 
of a new way of doing pQCD scattering am- 
plitude calculations that could have a dra- 
matic impact on QCD phenomenology. I will 
also discuss some important and novel LHC 
processes where more calculational effort is 
needed. Lack of space prevents a detailed 
review of all the interesting experimental re- 
sults on QCD-related processes presented at 
this Conference, but many of these are cov- 
ered in other plenary talks. 

2 Status of pQCD calculations 

For a broad class of 'hard' high-energy pro- 
cesses involving hadrons (i.e. suitably inclu- 
sive, with at least one large momentum trans- 
fer scale), the cross section can be calculated 
in QCD perturbation theory: 1 
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The classic example is the three-jet cross sec- 
tion in e + e~ — > hadrons, for which N = 1. 
The A, C\ and C2 coefficients define the LO, 
NLO and NNLO perturbative contributions 
respectively. Where the hard process involves 
two large but unequal energy scales, the per- 
turbative coefficients may be dominated by 




Figure 1. Interfacing a NLO pQCD calculation with 
parton showering in the production of a W or Z boson 
in hadron-hadron collisions. 

large logarithms of the two scales, 

da = Actg [ 1 + (ci,iZ + Ci >0 )q!s 

+ (c 2 , 2 i 2 + c 2 .iL + +c 2:0 )a| + •■• ] (2) 

where L = log(M/gr), log(l/x), log(l — 
T), ... 3> 1, for example. In this case one can 
often resum the leading and some subleading 
logarithms to all orders to improve the per- 
turbative prediction, the c n<n , c nj „_i, c n , n -2 
coefficients defining the LL, NLL, NNLL con- 
tributions respectively. We will show exam- 
ples of state-of-the-art resummed pQCD phe- 
nomenology below. 

Perturbative QCD calculations become 
more complex the more external quarks 
and gluons are involved. Thanks to the 
development of automated codes such as 
MADGRAPH 2 and CompHEP 3 , in princi- 
ple any tree-level (leading-order) scattering 
amplitude can be calculated. However the 
large renormalization scale dependence gen- 
erated by the overall as{/J?) N factor means 
that precision predictions for overall rates 
are not possible at this order. This is par- 
tially solved by extending the calculation to 
NLO, where the coefficient C\ receives real 
and virtual contributions: d<jy^ + da^ +1 \ 
Although scale dependence is reduced (for a 
cross section calculated in pQCD at order 
a%, ii 2 da N /d(i 2 = 0(a* +1 )) it still domi- 
nates the uncertainty on many NLO as mea- 
surements. Another improvement at NLO is 
that jet structure begins to emerge, i.e. a 
jet can contain one or two partons and there- 
fore has a perturbative 'width', in contrast to 



LO where jet = parton. Indeed an important 
advance in NLO technology has been the de- 
velopment of techniques for interfacing the 
perturbative NLO calculation with parton 
shower models (i.e. leading logarithm par- 
ton branching followed by hadronization), il- 
lustrated schematically in Fig. 1. The proce- 
dure for doing this is highly non-trivial — one 
must avoid double counting gluon emissions, 
and deal with cancelling negative and posi- 
tive singularities from the two types of fixed- 
order contribution shown in { } in Fig. 1. 

The most advanced formalism of this 
type is the MC@NLO programme of Frix- 
ione et a/. 4 , which combines the full pQCD 
NLO calculation with the HERWIG parton 
shower Monte Carlo. MC@NLO already 
contains a number of important processes 
for hadron collider phenomenology, includ- 
ing W, Z, WW, WZ, ZZ, bb, tt and H 
production. The benefits of combing fixed 
order N"LO and parton showers are obvi- 
ous: one retains the correct overall rate 
(with reduced scale dependence) and the full 
hard-scattering kinematics, while generating 
a complete event picture and a consistent 
treatment of collinear logarithms to all or- 
ders. This is illustrated in Fig. 2, which shows 
the MC@NLO prediction for the pt distribu- 
tion of tt pairs at the Tevatron pp collider. 
One would expect the fixed-order and the 
parton shower (resummed leading logarithm) 
approaches to give the correct prediction at 
low and high px respectively, and indeed the 
MC@NLO prediction interpolates smoothly 
between them. 

In contrast to LO, there is as yet no au- 
tomated technique for NLO calculations - 
each process has to be treated individually 
- although this is currently a topic of de- 
tailed study. As a result, there are many 
processes for which the full NLO corrections 
are not yet known. The most phcnomcnolog- 
ically important are those involving the pro- 
duction of multiple gauge bosons and heavy 
quarks in hadron collisions, which form the 
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Figure 2. MC@NL0 4 prediction for the tt transverse 
momentum distribution at the Tevatron collider. 



background to many New Physics processes. 5 
For example, the 0(a s ) qq,gg — » ttbb pro- 
cess forms an irreducible background to as- 
sociated tiH(—> bb) production at the LHC. 
The signal to background ratio is 0(1), but 
the latter is currently known only at LO and 
there is therefore a large scale dependence in 
the prediction. 
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Figure 3. Summary of as(M^) measurements, from 
S. Bethke. 6 



3 as measurements 

Figure 3, from S. Bethke, 6 summarizes the 
as(Mz) measurements from some of the 
most accurate recent determinations. The ta- 
ble contains a mix of NLO and NNLO mea- 
surements. For experiments performed at en- 
ergy scales different from Mz, the as val- 
ues measured at fi 2 — Ql xp are converted 
to as(M%) using the standard expressions. 1 
The consistency of the various measurements 
is remarkable - - a s is indeed a univer- 
sal parameter! Defining a 'world average' 
value presents a technical difficulty, however. 
Since the errors on most of the measure- 
ments are largely theoretical — often based 
on estimates of unknown higher-order correc- 
tions or non-perturbative effects — and nei- 
ther Gaussian nor completely independent, 
the overall error on the combined value of 
as(M§) cannot be obtained from standard 
statistical techniques, see Ref. 6 for details. 



Bethkc's (2004) world average value and er- 
ror is 

a MS,NNLO( M 2) = Q.1182 ± 0.0027 (3) 

In view of the consistency of all the measure- 
ments, and in particular of those with the 
smallest uncertainties, it seems unlikely that 
future 'world average' values of as will de- 
viate significantly, if at all, from the current 
value given in (3). Indeed, the correspond- 
ing world average value in 2002 was almost 
identical at 0.1183 ± 0.0027. 6 

At this Conference, a number of new 
as measurements have been reported, from 
processes including (i) HERA jet cross sec- 
tions and shape variables, (ii) 4-jet shape and 
jet shape moments from reanalysed JADE 
(e + e~) data, and (iii) LEP jet shape observ- 
ables and 4-jet rates. All these are NLO mea- 
surements, and all are consistent with the val- 
ues shown in Fig. 3. 
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Table 1. Summary of quantities known (exactly) to NNLO in pQCD. References are given to recent calcula- 
tions only. Not included are other partial or approximate (e.g. soft, collinear) results and NNLL improvements. 



4 NNLO 

As experimental present and future measure- 
ments reach the few per cent accuracy level, 
pQCD calculations at NNLO or higher are 
required. NNLO is in fact the current pertur- 
bative frontier, although NNLO corrections 
for a number of processes, particularly inclu- 
sive quantities with an external electroweak 
gauge boson providing the large energy scale, 
for example the e + e~ — » hadrons total cross 
section, have been known for some time. 1 

Table f summarises the quantities for 
which the NNLO pQCD corrections are cur- 
rently known. The calculations of those pro- 
cesses with references listed have all been 
completed within the past two years. Per- 
haps the most important NNLO calcula- 
tion still outstanding is the inclusive high- 
Et jet distribution at hadron colliders. This 
is needed to complete the ingredients for a 
full NNLO parton distribution function (pdf) 
global analysis (see below) , and to search for 
a New Physics contribution to the high-Er 
tail of the distribution. Schematically, 

J — jet 

— - = a A + a 3 (B + 2b LA) + a 4 (C 
dhjT 
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Figure 4. Dependence of the predictions for the Teva- 
tron jet inclusive cross section at Et = 100 GeV on 
the renormalization/factorization scale fi in the MS 
scheme, from Ref. 21 

+3b LB + (3blL 2 + 2b 1 L)A) (4) 

with a = for renormalization and 

factorization scale Although the NNLO 
correction coefficient C is not yet known, 
its likely effect in reducing factorization and 
renormalization scheme dependence is illus- 
trated in Fig. 4 21 , where the cross section is 
plotted for the choices C = 0, ±B 2 /A. 

The complexity of the NNLO jet calcu- 
lations stems from the fact that the singu- 



%(x,Q 2 ) 



»W -Si J {^M^Q*) + P..(*«.to<f.«")} 

dg(x,Q 2 ) a s f 1 dy r x 2 a; 2 "i 



(5) 



lar parts of the different components of the 
0(a s ) calculation (the 2-loop, 2-parton fi- 
nal state; the 1-loop-squared, 2-parton fi- 
nal state, the 1-loop, 3-parton final state; 
the tree-level, 4-parton final state) have to 
be identified and cancelled, leaving behind 
a non-singular net contribution that can be 
evaluated numerically to yield the coefficient 
C in (4). The key is to identify and calcu- 
late the various 'subtraction terms' that add 
and subtract to render the loop (analytically) 
and real emission (numerically) contributions 
separately finite, see the review by Glover 21 
for a more detailed discussion and list of ref- 
erences to recent work in this area. 

5 Three-loop splitting functions 

Hard-scattering cross sections in hadron- 
hadron collisions (e.g. pp — > H + X) are 
obtained by convoluting subprocess cross sec- 
tions (a) with parton distribution functions, 
fi(x,Q 2 ), whose factorization-scale (Q 2 ) de- 
pendence is determined by the DGLAP evo- 
lution equations (5). Consistency requires 
the subprocess cross sections and DGLAP 
splitting functions to be calculated to the 
same order in perturbation theory: 

P(x, a s ) = P (0) + a s P {1) + a|P< 2 > + ... 

a = a^ +as& W +a 2 a^ + ... ( 6 ) 

The full set of LO and NLO splitting func- 
tions were calculated in the 1970s and 1980s 
- see Ref. 1 for details. Calculations of the 
NNLO coefficient functions of the various 
hadron collider subprocess cross sections date 
back to the early 1990s, see Table 1. However 
it was only a few months ago that the full set 
of NNLO ('three-loop') splitting functions fi- 
nally became available when, in a landmark 
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Figure 5. One of the three-loop (NNLO) non-singlet 
splitting functions from Moch, Vermaseren and 
Vogt. 7 ' 8 . Sec text for details of the various curves. 



calculation, Moch, Vermaseren and Vogt 7,8 
(MVV) completed the determination of all 
the components of the full splitting function 
matrix at NNLO, i.e. the pjj%, P$, P^} 
and P^ functions of (5), (6). Both (very 
lengthy) complete analytic forms and (com- 
pact) numerical approximations for practical 
applications are provided. One of the three 
-loop (NNLO) non-singlet splitting functions 
from MVV is shown in Fig. 5. In the upper 
figure, the exact function (solid line) is com- 
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Figure 7. Resummed pQCD prediction for the Higgs 
transverse momentum distribution at the LHC, from 
Bozzi et al. 25 



perturbation series is well under control, and 
there is also good agreement between theory 
and experiment. 



Figure 6. Predictions for the W and Z total cross 
sections at the Tevatron and LHC, from MRST. 23,24 



pared with previous approximate forms 22 fit- 
ted to a finite number of exact moments and 
known x — > 0, 1 leading behaviour. In the 
lower figure, the small- a; leading-logarithm, 
next-to-leading-logarithm, etc. approxima- 
tions are compared with the exact result. 

Using the results of MVV, consistent 
NNLO calculations of hard-scattering cross 
sections at hadron colliders can now be per- 
formed. One of the most important calcu- 
lations is the prediction for the total W or 
Z production cross sections. With sufficient 
theoretical precision, these could be used as a 
luminosity monitor at the Tevatron or LHC. 
A recent NNLO calculation from Martin et 
al. (MRST), 23,24 which uses the full MVV 
three-loop splitting functions, is shown in 
Fig. 6, together with recent CDF and DO 
Tevatron cross-section measurements. Tak- 
ing all sources of uncertainties in the theo- 
retical analysis into account, MRST quote a 
±4% uncertainty in the prediction, 24 shown 
as the horizontal dashed lines bracketing the 
NNLO cross sections in Fig. 6. Evidently the 



6 Resummation 

Effort continues to refine the predictions for 
'Sudakov' processes, of which the most fa- 
miliar examples are the thrust distribution 
near T = 1 in e + e _ collisions, and the W, 
Z or Higgs transverse momentum distribu- 
tion at small qT in hadron-hadron collisions. 
In the latter case, a precise knowledge of 
the qt distribution is necessary in order to 
be able to accurately measure the proper- 
ties (mass, width, etc.) of the heavy bo- 
son. At small qr the perturbation scries 
becomes dominated by terms of the form 

E„, m < 2 „-i a s lo g m ( M 79T) and resumma- 
tion is necessary. To obtain a prediction valid 
over the full qT range, the resummed cross 
section must be matched onto the exact fixed- 
order calculation at large qT ■ This is not triv- 
ial, since the resummed calculation ignores 
power-law contributions of order (q T /M 2 ) N 
that contribute at large qT- Matching is 
achieved by first subtracting from the full 
fixed-order calculation those logarithmic con- 
tributions that form part of the resummed 
calculation, and then adding the two types of 
contribution. Figure 7 shows a recent calcu- 
lation of the Higgs qT distribution at the LHC 



+ + + + 

Figure 8. Vanishing n-gluon helicity amplitudes. 



from Bozzi et al. which combines the fixed- 
order (NLO) and resummed (NNLL) contri- 
butions. The effect of resummation is clearly 
seen at small qx in the left-hand figure. The 
right-hand figure shows that the scale inde- 
pendence of the fixed-order result is main- 
tained in the full calculation over the whole 
qx range. 

7 The CSW technique for QCD 
amplitudes 

Although the numerical calculation of tree- 
level QCD scattering amplitudes can be auto- 
mated, the method is very much 'brute force' 
and the complexity associated with multipar- 
ton amplitudes soon saturates the computer 
capability. And as already mentioned, equiv- 
alent automation for loop amplitudes seems 
some way off. Analytic expressions are of 
course only feasible for few-particle ampli- 
tudes, and lengthy analytic expressions are in 
any case of little practical use for numerical 
computations: see, for example, the expres- 
sions for the three-loop splitting functions in 
Rcfs. 7 < 8 . 

In the past year, a new technique based 
on an idea by Witten 27 and subsequently de- 
veloped by Cachazo, Svrcek and Witten 28 
(CSW) may be the long-awaited break- 
through in making QCD multiparticle tree- 
level and loop amplitudes tractable. The 
discussion of the CSW method starts with 
an observation by Parke and Taylor 29 in 1986 
concerning the helicity dependence of n-gluon 



"The Parallel Session contribution of Zvi Bern gives 
a very clear and comprehensive introduction to the 
CSW method. 
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Figure 9. The MHV amplitude. 



QCD tree-level scattering amplitudes. If all 
the external gluon hclicities are the same, or 
if only one has a different sign, then the am- 
plitude vanishes, sec Fig. 8. Although other 
amplitudes are non-zero in general, the Maxi- 
mum Helicity Violating (MHV) amplitude, in 
which two and only two external gluons have 
the opposite helicity (see Fig. 9) appears to 
have a special status in that it can be writ- 
ten in an extremely compact form in terms 
of spinor products: 



A 



MHV — I 9 S 



(r,. 



n-=iO-,i+i) 



(7) 



where colour factors have been suppressed, 
and the external momenta arc labelled i = 
1, ...,n with r and s the two lines with oppo- 
site helicity. No other non-zero helicity am- 
plitudes can be written in such a simple form. 
The ( ) quantities in (7) are spinor products 
defined by 



= u-(pi)u+(pj) , \ = v / 2p l ■ Pj 

(8) 

The spinor product is a complex func- 

tion of the two four-momenta pi and pj , the 
exact form depending on the choice of spinor 
representation. For example, a convenient 
choice for programming purposes is 30 
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for momentum 4- vectors p M = (p° ,p x ,p v ,p z ). 




Figure 10. Construction of a non-MHV amplitude 
from two MHV amplitudes and a scalar propagator. 

In the CSW method, the MHV ampli- 
tudes are elevated to effective vertices in a 
new scalar graph approach to tree-level am- 
plitude computation. Other helicity am- 
plitudes are obtained by stitching together 
MHV amplitudes with scalar propagators, 
taking account of all possible permutations 
of helicities. As an example, Fig. 10 shows 
how the ( h + + ...) amplitude is ob- 
tained by combining two ( h + +...) MHV 

amplitudes with a 1/q 2 scalar propagator. 
The result is a compact expression written 
in terms of spinor products (generalised to al- 
low for off-shell momenta) much simpler than 
the corresponding output obtained using the 
automated programmes. It is of course triv- 
ial to check that the CSW-method ampli- 
tudes agree with those obtained using stan- 
dard techniques. It is also possible to develop 
a recursive technique that allows MHV ver- 
tices to be stitched together in a systematic 
way in order to build up more complicated 
amplitudes. 31 

Two immediate questions are: (i) why 
does the method 'work', and (ii) can it be 
extended to more general QCD amplitudes, 
including those with loops? The answer to (i) 
lies in the fact that the perturbative expan- 
sion of J\f = 4 Supersymmetric Yang-Mills 
(SYM) gauge theory is equivalent to the in- 
stanton expansion of a certain type of string 
theory in supersymmetric twistor space. 27 . 
Tree- level amplitudes with n + 1 negative he- 
licity gluons are related to D— instantons of 
degree n. The MHV vertices are localised 
on a line in twistor space, and the lines are 



connected by off-shell propagators. Because 
of this topological structure, amplitudes in 
unbroken gauge theories are therefore much 
simpler than their Feynman diagram expan- 
sion would suggest. Since N = 4 SYM gauge 
theory shares the same gauge boson ampli- 
tudes as pure Yang-Mills theory, the results 
for gluonic amplitudes in QCD follow imme- 
diately. 

Since the original CSW paper, there has 
been a steadily increasing amount of theoret- 
ical ac tivity. 31 < 32 ' 33 < 34 < 35 < 36 < 37 The method has 
been extended to amplitudes with fermions 32 , 
and the important question of whether there 
is any equivalent simplification for QCD am- 
plitudes containing loops has begun to be 
addressed 36 ' 37 . While the MHV method ap- 
pears to work at the one-loop level in Af = 4 
SYM, the extension to loop amplitudes in 
non- supersymmetric theories auch as QCD 
appears to be more problematic. Neverthe- 
less given the current level of theoretical ac- 
tivity, there is likely to be much progress in 
the near future. 

8 Diffractive Higgs production and 
related processes 

The discovery of the Higgs boson is obvi- 
ously one of the main goals of the LHC. For 
the 'standard' search scenario, pp —> HX, 
via gg — > H, qq — > WH, gg — > ttH etc., 
the theory rests on the solid foundation of 
the QCD factorization theorem, and the to- 
tal rates and kinematic distributions can be 
accurately predicted for a given Mh- There 
is, however, no single optimum detection pro- 
cess but rather a range of possibilities de- 
pending on Mh, none of which is compelling 
on its own. It is important, therefore, to 
explore other 'non-standard' search scenar- 
ios. One that has received much attention 
recently (for a review of recent work and a 
list of references, see Ref. 38 ) is the exclusive 
process pp — > p © H © p, where the © sig- 
nals the presence of a rapidity gap. Provided 



that appropriate forward proton taggers can 
be installed, it may be possible to measure 
the Mh to high precision via the missing 
mass associated with the forward protons. 
This precision would then enhance the sig- 
nal over the continuum background, for ex- 
ample in the H — > bb channel. Further- 
more, the signal to background is also en- 
hanced by a J 2 = 0, P-even selection rule 
that suppresses the gg — > bb background. 39 ' 40 
Although the Standard Model Higgs boson 
is the obvious candidate for a new particle 
search in this channel, other exotic produc- 
tion processes are possible. Examples are 
SUSY Higgs bosons, gluino bound states, 
gravitons etc., indeed anything that couples 
strongly to gluons and survives the selection 
rule. 41 

However the calculation of processes such 
as pp — ► p®H(Bp presents a real challenge for 
the theory, since it involves both perturbative 
and non-perturbative aspects. In particular, 
one needs techniques for calculating not only 
the production amplitude shown in Fig. 11, 
but also the rapidity gap 'survival probabil- 
ity', i.e. the probability that strong inter- 
actions between the upper and lower proton 
systems in Fig. 11 do not populate the ra- 
pidity interval between them with additional 
soft hadrons. Preliminary indications are 
that signal to background ratios greater than 
one can indeed be achieved, 41 and the pos- 
sibility of measuring similar processes at the 
Tevatron (for example, pp — > p © \c ® P) of- 
fers a means of checking and calibrating the 
Higgs calculation. On the experimental side, 
the missing mass resolution is crucial with 
AM m ; ss ~ 1 GcV being the goal, see for ex- 
ample the recent study of Rcf. 42 . 
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